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ABSTRACT 

High resolution spectral observations of the polar EF Eridani obtained in a low state with X-Shooter revealed narrow emission lines from 
the irradiated secondary. The lines were most prominent in the near-lR Can-triplet, the more intensive Ha line had additional emission likely 
originating from an accretion stream. The lines with a radial velocity amplitude, K'2 = 385 ± 4 km s"', serve as tracer of the otherwise unseen 
companion. The mass function implies a massive white dwarf with M„j> 0.65 M© at 3cr confidence, and a short distance to the binary, 
(i ~ 111 pc (< 145 pc at 3iT confidence). The spectral energy distribution from the UV to the IR together with the high mass ratio gives further 
strong evidence of EFEri being a post period-minimum object with M2 < 0.06 Mq. 
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1. Introduction 

EF Eridani was only the third Polar (AM Herculis binary) dis- 
covered after AM Herculis itself and VV Puppis and the first 
one discovered initially as a bright X-ray source with spec- 
troscopic and polarimetric identification afterwards (Williams 
et al. 1979; Griflaths et al. 1979; Tapia 1979). The long-term 
light curve assembled from the Harvard plate stacks displays 
variability of the binary between Z? = 15 - 18 mag (Griffiths 
et al. 1979), naturally explained as high and low accretion 
states. Since its discovery in 1 978 it remained in a high state un- 
till 1997 (Wheatley & Ramsay 1998); the current low state was 
interrupted by just three short-lasting accretion events (Mukai 
et al. 2008). 

During high and low states EFEri behaved as the text- 
book example of AM Herculis binaries (Beuermann et al. 1987; 
Schwope et al. 2007), at least as far as energy conversion and 
release in the accretion column was concerned. It changed from 
shock-heating/bremsstrahlung cooling in the high state to parti- 
cle heating/cyclotron cooling in the low state, when the overall 
accretion rate is reduced by a factor 100 or more. 

The high state emission line spectrum of EF Eri was stud- 
ied by Verbunt et al. (1980); Bailey & Ward (1981); Hutchings 
et al. (1982); Mukai & Charles (1985). They describe a com- 
plex multi-component line pattern of high and low ionization 
species, nevertheless stable from epoch to epoch. None of the 

* Based on observations made with ESO Telescopes at the Paranal 
Observatory under programme ID 60.A-9436(A) 



several line components could be uniquely associated with a 
physical structure in the binary apart from the high-velocity 
fine wings which clearly originate in streaming matter between 
the two stars. 

The low state gave view to the component stars in 
EFEri. A complex magnetic field topology for the white 
dwarf with strong quadrupolar components was infered from 
Zeeman tomography based on low-resolution VET/FORS spec- 
tra (Beuermann et al. 2007). Evidence for a high-mass white 
dwarf, M„d= 0.82- 1.00 Mq was derived by fitting the GINGA 
X-ray spectrum (Cropper et al. 1998). Since the model of 
Cropper et al. in some cases and for unknown reasons revealed 
masses in stark contrast to mass estimates based on radial ve- 
locities or spectroscopic parallaxes (see e.g. Gansicke et al. 
1995; Catalan et al. 1999, for AM Her and QQ Vul, respec- 
tively) independent confirmation is necessary. Recently Howell 
et al. (2006) presented spectroscopic monitoring of EFEri with 
the 1 .3m SMARTS telescope and they report the detection of 
several lines in the low state which they ascribe to chromo- 
spheric activity of the secondary. Ha was persistently observed 
throughout the orbital cycle. A sine fit to the Ha radial veloc- 
ity curve revealed K2 = 269 + 18 km s~\ regarded by them as 
projected orbital velocity of the secondary. 

EFEri is a key object for binary stellar evolution with 
high-field white dwarfs. Close to the CV period minimum, 
f orb = 81 min, it may harbour a late-type main-sequence star 
or, as a period bouncing system, a brown dwarf-like degenerate 
object as secondary star. Evidence for a sub-stellar secondary 
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Fig. 1. Trailed spectrograms of EFEri of tlie Ha (left) and Can 8542 emission lines (right) obtained August 13, 2009, in a low state. The 
additional line in the right panel is Can 8498. 



was derived from an upper limit TiO band strength in the first 
low-state NTT spectra (Beuermann et al. 2000). Likely spec- 
troscopic and photometric signatures were reported by Howell 
& Ciardi (2001) and Harrison et al. (2003) but the claims of the 
detected secondary remained unconfirmed. Infrared cyclotron 
emission from a tenuous accretion plasma is definitely con- 
tributing to the excess infrared emission, the contribution of 
the low mass donor is still uncertain (Harrison et al. 2004; 
Schwope et al. 2007; Campbell et al. 2008). 

To search for lines useful for a radial velocity study and 
establish the spectral energy distribution from the ultraviolet 
to the infrared we proposed to include EFEri in the Science 
Verification programme for the VLT 2nd generation instru- 
ment X-Shooter (D'Odorico et al. 2006). The observations 
took place on August 13, 2009, and the raw data were made 
available on August 25. We take the opportunity to single out 
one important result from an early pipeline reduction of the 
complex data set, highlighting the capability of the instrument 
for binary star research and addressing an important issue in 
EFEri, its radial velocity solution, and the mass constraints 
which rest on it. 



2. X-Shooter spectroscopy of EF Eri 

EF Eri was selected for X-Shooter Science Verification obser- 
vations for a total of 1*. One observation block was created to 
be executed in service mode. Telescope and instrument over- 
heads reduced the time on target to 2080 s, which was split in 
equal parts to obtain 1 1 spectral exposures of length 189 s each. 
The observations started August 13, 2009, at UT 07:45:32, and 
covered binary cycle 197488.508 - 197489.076 using the up- 
dated ephemeris of the X-ray dip Beuermann et al. (2007) and 
the period derived by Piirola et al. (1987). A phase resolution 
of 0.056 of the f orb = 4861 s binary was achieved. 

The seeing was slightly variable between 0'.'83 and r.'02 
thus always almost matching the chosen slit width of V.'2. 
Spectrum #7 at phase 0.46 was obtained under the least op- 
timum seeing, the continuum and line flux suffer slightly from 
likely slit loss (see Figs. 1 and 3). 



The X-shooter data were reduced with a pre-release version 
of the pipeline (Goldoni et al. 2006). Data from the three arms 
were reduced with similar procedures. The pipeline uses a stan- 
dard bias subtraction in the UVB and VIS arms, and dark frame 
subtraction in the NIR arm. Sky emission lines are subtracted 
from the curved orders using a resampling procedure (Kelson 
2003). The orders are traced, extracted and wavelength cali- 
brated and rectified. The final results provide frames of both 
individually extracted orders and their associated errors, and 
a frame with all the orders merged weighted according to the 
errors. 

One-dimensional spectra were finally extracted with stan- 
dard MIDAS procedures. The spectra were binned on a 0.25 A 
grid in the visual arm, the spectral resolution was estimated 
from the FWHM width of night sky Unes to R = 6000 
at 8500 A, typical uncertainties of the wavelength calibration 
were determined to be 0.1 - 0.2 pixels. No science-grade spec- 
tral response function for our observing night was available 
at the time of writing. However, a quantitative analysis of the 
main etnission lines is nevertheless possible currently and we 
concentrate in this paper on the spectral region around Ha and 
the Can-triplett which are covered by spectral orders 25 and 19 
of X-Shooter's visual arm, respectively. 

Continuum flux variability over the 45 min of observations 
was found to be below 7% at short wavelengths, < 7200 A. 
This low level could already be caused by slit losses due to 
variable seeing and is not necessarily intrinsic to the source. A 
red spectral component longward of ~7200 A likely similar to 
that observed by Wheatley & Ramsay (1998) was present in 
the first 2 spectra. Pronounced Zeeman absorption lines from 
the white dwarf and the low level of orbital variability clearly 
reveal that EFEri was observed in a low state of accretion. 

The 1 1 spectra were arranged in a two-dimensional image 
to mimic a trailed spectrogram. A few emission lines could 
be recognized, most prominently Ha, Hel5875 and the Can- 
triplet (8498, 8542, 8662). An average spectrum around the 
Can-triplet is shown in Fig. 2 and cutouts of the trailed spec- 
trogram around Ha and the most prominent Can line at 8542 A 
are displayed in Fig. 1. 
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Fig. 2. Mean X-Shooter spectrum of EFEri (visual arm, 19th order) of 
exposures #3 -#9. The individual spectra were transformed to zero ra- 
dial velocity before averaging. Observed Can and possible H-Paschen 
lines are indicated. 
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Fig. 3. Radial velocity and flux variations of the Can 8542 emission 
line. The dashed line in the upper panel indicates the best fitting 
sine curve and dashed lines in the lower panel give predicted line 
flux variations based on our irradiation model for orbital inclinations 
i - 20°, 30° . . . 90° (from top to bottom). 



While Ha has a complex line profile with a sharp compo- 
nent varying sinusoidally and line wings extending to several 
hundred kilometers away from the naiTow core, the Can lines 
display just the naiTow line component. It displays pronounced 
flux variability with maximum flux occurring around orbital 
phase 0.5 (phase convention derived below). Giving the striking 
similarity to line profile variability of other polars observed in 
high and low states (see the collection of Doppler maps in e.g. 
Schwope 2001, and references therein) the narrow line compo- 
nent can unequivocally be assigned to the secondary star. The 
Ha line has additional components, perhaps best visible in the 
first spectrum at phase 0.12. Another rather narrow component 
offset from the sine curve of the naiTow emisison line can be 
recognized therein. 

We measured radial velocities of the Can 8542 line by fit- 
ting single Gaussians to the observed line profiles. We focus 
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Fig. 4. Radial velocity factor K2lK'^ to correct from center of light 
to center of mass based on our irradiation model for a star filling its 
Roche lobe. 



here on the brightest of the three lines for signal/noise con- 
siderations. A line position could not be measured in the first 
spectrum due to the faintness of the line. In Fig. 2 the positions 
of the Can lines and of the Hydrogen Paschen lines are indi- 
cated. The radial velocity corrected average spectrum displays 
weak Paschen emission lines of P9, PI 1, and P14, respectively. 
Hence, the blending lines PIO, PI 2, and PI 3 are weakly con- 
tributing to the Can lines but due to their faintness cannot be 
recognized in the individual spectra and are not considered in 
our analysis. 

Results of radial velocity measurements of the Caii8542 
line are shown in Fig. 3. An unweighted sine fit to the ra- 
dial velocity curve revealed a radial velocity amplitude K'2 - 
385 + 4 km s ' . The epoch of blue-to-red zero crossing at BJD 
= 2455056'?8188 ± 0!0003 determines orbital phase = 0.0 
throughout this paper and is a robust measurement of infe- 
rior conjunction of the secondary in EFEri. Deviations from 
the sine curve are negligible. The line flux is highly variable, 
reaching a maximum plateau between phases 0.4 and 0.55 
and declines before and thereafter. Similarly variable is the 
FWHM line width which is ~ 130 km s ' in the phase interval 
(f) - 0.4 - 0.65 (i.e. the lines are resolved in this phase interval) 
and gradually declines down to the instrumental resolution of 
~50km s ' observed in spectrum #2 at phase 0.17. 

2.1. Results and analysis 

The line profile variations of the Can lines leave little doubt 
about their origin. They were emitted from the hemisphere of 
the secondary facing the white dwarf and being iiTadiated by X- 
ray and UV radiation from the accretion spot and white dwarf. 
This will be our working hypothesis in the following. 

To constrain the mass of the white dwarf we need to es- 
timate the true orbital velocity (corrected for iiTadiation and 
inclination) and the mass of the secondary. 

We use an irradiation model (Beuermann & Thomas 1990) 
to coiTect the observed K2 from center-of-light to center-of- 
mass and by that predict the true orbital velocity. The result- 
ing correction factor is a function of the mass ratio in the first 
place. The very weak dependance on the orbital inclination is 
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Fig. 5. (a) White dwarf mass function based on our irradiation model. 
Solid lines are for an 0.045 M©, dotted lines for an 0.072 M© sec- 
ondary. The three lines each illustrate the 3(t range around the nominal 
velocity K'^ = 385 km s"'. (b) White dwarf distance in EFEri for as- 
sumed temperatures of 10000 K (black solid line) and 1 1000],K (green 
dashed line) (c) Maximum J-band absolute brightness of the secondary 
in EFEri for the given distance. The black solid line was computed for 
a J-band flux from the secondary of 0.09 x 10"'* erg cm"^ s"' A"', the 
dashed red line for 0. 1 8 x 1 0" ' * erg cm"^ s" ' A" ' . Within the shaded re- 
gion the transition from a stellar (left) to a sub-stellar (right) secondary 
(left) occurs according to Knigge (2007). 



negligible. The correction factor is shown in Fig. 4, it is for all 
acceptable mass ratios larger than 15%. 

The orbital inclination of EFEri was constrained by the 
occuiTence of X-ray absorption dips to / > S. The angle 6 
is defined by the local magnetic field vector in the accretion 
spot on the surface of the white dwarf (see e.g. Beuermann 
et al. 1987). The inclination was further constrained by cy- 
clotron models applied to photoelectric polarimetry (Cropper 
1995; Piirola et al. 1987), and by modeling X-ray and UV Hght 
curves (Beuermann et al. 1987; Schwope et al. 2007). All au- 
thors agree that the inclination is most probably between 50° 
and 65°, although values as high as / = 75° cannot be com- 
pletely excluded. 

The inclination influences the visibility of the irradiated re- 
gion of the secondary and thus gives an imprint on the line flux 
variability. It can be estimated from the present observations. 
We use the iiTadiation model to calculate Doppler-broadened 
line profiles and generate emission-line light curves as a func- 
tion of the orbital inclination. Results are shown in Fig. 3 with 
a normalization chosen by eye, the same value for all light 
curves. Line flux variability rules out / < 40° in agreement with 
previous estimates. The phase coverage and the signal-to-noise 
ratio of the present data is insufficient to draw more stringent 
constraints on /. There are differences between the model and 
the data around maximum phase, e.g. the models do not show 
the plateau between <p = 0.4 - 0.6. The slight depression of the 
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Fig. 6. Spectral energy distribution of EFEri. For an explanation of 
the used symbols and the data displayed see text. 



line flux in spectrum #7 at phase 0.46 is likely due to slit losses 
and thus apparent only. 

The mass of the secondary could not too strongly con- 
strained observationally in the past, M2 < 0.142 (see the 
discussion in Beuermann et al. 2000). We thus make use of 
Knigge's donor sequence (Knigge 2006, 2007) to further con- 
strain the likely mass of the white dwarf. At the orbital period 
of 4861 s close to the CV period minimum the object might be 
a pre-period minimum ordinary CV or a post period minimum 
object, a so-called period bouncer. Knigge's relation predicts 
M2 = 0.072 Mo for the first, M2 = 0.045 Mq for the second 
case. Using the irradiation model we calculated as a func- 
tion of the mass ratio and of the inclination for the two assumed 
masses. The curves in Fig. 5(a) imply Mwd> 0.66 M©, unless 
the secondary is M2 < 0.045 Mq. 

Thorstensen (2003) calculated the distance to EFEri by 
combining information from parallax observations with a num- 
ber of other pieces of information, including an estimate of 
the a-priori velocity distribution towards EF Eri. The constraint 
from parallax for EF Eri was weak, so the priors dominated his 
Bayesian calculation. This lead Thorstensen to give two dis- 
tances, one with a velocity prior including a high-velocity pop- 
ulation, leading to a large distance of 163 pc, and one without 
this population, leading to a short distance of 113 pc. 

Inference on the likely distance can be made from the ob- 
served brightness and the derived temperature of the white 
dwarf. The effective temperature of the white dwarf is likely 
below lOOOOK (Schwope et al. 2007, Szkody et al. 2010 sub- 
mitted), although a temperature as high as llOOOK is not 
completely ruled out. The implied distances for the two as- 
sumed temperatures are illustrated in Fig. 5(b). A robust up- 
per limit for the distance to the white dwarf assuming 1 lOOOK 
is 168 pc, at lOOOOK the upper limit distance is 144 pc, for 
an inclination of 60° the distance is 1 1 1 pc and 96 pc for 
degenerate/non-degenerate secondaries, respectively (assum- 
ing Teff = lOOOOK). Fig. 5(a) also reveals that the orbital incli- 
nation has to be larger than 44° to not exceed Chandrasekhar's 
mass limit. 
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Taking M„d= 0.9 Mg at face value the inclination becomes 

60° or 66°(for a degenerate/non-degenerate secondary), in ex- 
cellent agreement with optical/UV light curve modeling. 

In Fig. 6 we show an updated version of the ultra vio- 
let/optical/IR spectral energy distribution of the binary (up- 
date with respect to Schwope et al. 2007). Apart from an op- 
tical low-state spectrum adapted from Harrison et al. (2004) it 
comprises GALEX photometry at orbital minimum and maxi- 
mum (Szkody et al. 2006), XMM-OM optical and visual pho- 
tometry (Schwope et al. 2007, UV red dots), orbital phase- 
resolved infra-red photometry in JHK (Harrison et al. 2003) 
and SPITZER infrared data from Brinkworth et al. (2007) (for 
a Spitzer spectrum and a further discussion of the mid-lR 
spectral energy distribution see Hoard et al. 2007, data points 
shown as blue dots). The green lines illustrate the combined 
white-dwarf/accretion spot model adapted from Schwope et al. 
(2007). 

Extrapolation of the white-dwarf/accretion-spot model into 
the infrared spectral regime shows that the most sensitive limits 
on the secondary can likely be derived from the J-band. The H 
and K bands have a much stronger dust component (Hoard et al. 

2007) . 

The low-state infrared spectrum and variability is dom- 
inated by cyclotron emission at all phases (Campbell et al. 

2008) . The observed J-band minimum flux corrected for 
the contribution from the white dwarf thus represents an 
upper limit to the flux of the secondary, = 9 x 
10"^^ergcm~-^s~^ A~^(ffiy = 18.87). We searched Knigge's 
donor sequence for a secondary that does not violate the bright- 
ness limit for a given distance. This limit, the maximum possi- 
ble absolute J-band magnitude of a donor, is shown in Fig. 5(c) 
as a function of the distance. In an erratum to the original publi- 
cation Knigge (2007) gives Mj = 13.2 at the stellar/sub- stellar 
boundary. However, at period minimum his IR-sequences show 
wiggles due to switches in the model atmosphere grids used. At 
the stellar/sub- stellar boundary Mj might be as bright as 12.6 
if one interpolates the tabulated values with a smooth func- 
tion. The region of uncertainty is indicate as a shaded box in 
Fig. 5(c). Fortunately, at P = 81 min we are away from this un- 
certain region on either side (stellar or sub-stellar). Hence from 
Fig. 5(c), the secondary has to be sub-steUar at any distance 
shorterthan~140(or~180pc,iftheboundaryisatMy = 12.6). 

According to Knigge's donor sequence, a stellar secondary 
at the orbital period of 81 min has Mj - 12.05, and the com- 
bined flux of the white dwarf (including the accretion spot) plus 
such a secondary would not exceed the observed J-band min- 
imum flux only at a distance of 230 pc (blue line in Fig. 6). 
Hence, if Knigge's sequence is applicable to EFEri and if the 
observed J-band minimum flux is a valid proxy of the sec- 
ondary's brightness, a stellar secondary in EFEri can likely be 
ruled out. 

A substellar Knigge-donor plus white dwarf will not ex- 
ceed mjimin) at a distance of 135 pc (magenta line in Fig. 6). 
This value is rather close to the mentioned boundary of 140pc 
and due to the non-monotonic behaviour of My around period 
minimum in Knigge's tables. 



3. Discussion and outloolc 

Using X-Shooter and exposing just 35 min on target we found 
good evidence for a massive white dwarf in EF Eri and a sub- 
stellar secondary, establishing the object as a period bouncing 
CV at a distance of about llOpc. The clues are based on line 
profile variations of the near-IR Can triplet which we found in 
accord with irradiation of the still unseen donor. 

We note that even our observed radial velocity amphtude, 
K^, which we regard as originating from the irradiated sec- 
ondary, is by 1 14 km s"^ larger than the Ha velocity pubUshed 
by Howell et al. (2006) and interpreted as true K2. This makes 
the discussion of the stellar masses presented there shaky. A 
Doppler tomographic analysis of the Can and the Ha fine of 
the X-Shooter data reveals the same velocities for the narrow 
emission line components. It appears likely that the Ha line 
observed by Howell et al. (2006) has contributions from other 
locations in the binary system, e.g. the stream. 

Our mass function was derived for two assumed secon- 
daries from Knigge's CV donor sequence. Investigating period- 
bouncing, eclipsing, non-magnetic CVs Littlefair et al. (2008) 
found deviations between those predicted and the observed 
masses. The secondaries were found to be more massive than 
predicted. The same may apply here. We thus underline the im- 
portance of an improved parallax for EF Eri to better constrain 
the mass, hence radius of the white dwarf as well as an im- 
proved determination of the orbital inclination (e.g. via Can 
line flux measurements with higher signal-to-noise and fuU- 
phase coverage) to better constrain the mass ratio. 

Nevertheless, the Ukely white dwarf mass found here of 
around 0.9 M© for / = 60° is in accord with X-ray spectral fits 
(Cropper et al. 1998). Such a high mass would conform with 
the mean mass of magnetic white dwarfs, 0.93 Mo(compared 
to 0.57 Mq for non-magnetic white dwarfs), as found by 
Wickramasinghe & Ferrario (2005). 

The direct discovery of the secondary in any short-period 
magnetic CV is an outstanding issue. Direct measurement or 
indirect evidence as presented here are very much important 
to compare the fraction of period-bouncing to pre-period min- 
imum objects and thus test population synthesis predictions 
(Howell et al. 2001) and to test the CV donor sequence at the 
bottom end. Further X-Shooting of carefully selected polar sys- 
tems is thus highly desirable. 
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